Search for: All records

Abstract. We present a newly developed upper-thermocline, open-ocean biogeochemical flux model that is complex and flexible enough to capture open-ocean ecosystem dynamics but reduced enough to incorporate into highly resolved numerical simulations and parameter optimization studies with limited additional computational cost. The model, which is derived from the full 56-state-variable Biogeochemical Flux Model (BFM56; Vichi et al., 2007), follows a biological and chemical functional group approach and allows for the development of critical non-Redfield nutrient ratios. Matter is expressed in units of carbon, nitrogen, and phosphate, following techniques used in more complex models. To reduce the overall computational cost and to focus on upper-thermocline, open-ocean, and non-iron-limited or non-silicate-limited conditions, the reduced model eliminates certain processes, such as benthic, silicate, and iron influences, and parameterizes others, such as the bacterial loop. The model explicitly tracks 17 state variables, divided into phytoplankton, zooplankton, dissolved organic matter, particulate organic matter, and nutrient groups. It is correspondingly called the Biogeochemical Flux Model 17 (BFM17). After describing BFM17, we couple it with the one-dimensional Princeton Ocean Model for validation using observational data from the Sargasso Sea. The results agree closely with observational data, giving correlations above 0.85, except for chlorophyll (0.63) and oxygen (0.37), as well as with corresponding results from BFM56, with correlations above 0.85, except for oxygen (0.56), including the ability to capture the subsurface chlorophyll maximum and bloom intensity. In comparison to previous models of similar size, BFM17 provides improved correlations between several model output fields and observational data, indicating that reproduction of in situ data can be achieved with a low number of variables, while maintaining the functional group approach. Notable additions to BFM17 over similar complexity models are the explicit tracking of dissolved oxygen, allowance for non-Redfield nutrient ratios, and both dissolved and particulate organic matter, all within the functional group framework. 

The distribution of oceanic biogeochemical tracers is fundamentally tied to physical dynamics at and below the mesoscale. Since global climate models rarely resolve those scales, turbulent transport is parameterized in terms of the large‐scale gradients in the mean tracer distribution and the physical fields. Here, we demonstrate that this form of the eddy flux is not necessarily appropriate for reactive tracers, such as nutrients and phytoplankton. In an idealized nutrient‐phytoplankton system, we show that the eddy flux of one tracer should depend on the gradients of itself and the other. For certain parameter regimes, incorporating cross‐diffusion can significantly improve the representation of both phytoplankton and nutrient eddy fluxes. We also show that the efficacy of eddy diffusion parameterizations requires timescale separation between the flow and reactions. This result has ramifications for parameterizing subgrid scale biogeochemistry in more complex ocean models since many biological processes have comparable timescales to submesoscale motions.